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- SUMMARY 



Vibrations of j?ofcary~wing aircraft may derive their 
energy from the rotation of tha rotor rather than from 
the air forces. A theoretical analysis of these vibra- 
tions is described and methods for its application are 
explained herein. 

The present paper also supersedes and extends the 
scope of the Advanot SU«tricte& Export entitled "Theory 
of Self-Excited Mechanical Oscillations of Hinged Rotor 
Elades," parts of which are in error. The theory has 
been extended to include the effects of unequal stiffness 
of the pylon for deflections in different directions and 
the effect of damping in the hinges and in the pylon. 
Both the derivation of the characteristic equation and 
the methods of application of the theory are given. In 
particular, the theory predicts the so-called "odd- 
frequency" self-excited speed range as well as the shaft- 
critical speed. Charts are presented from which the 
shaft-critical and the self-excited instabilities can be 
predicted for a great variety of cases. The influence of 
each physical parameter upon the instabilities has been 
obtained. The comprehensive treatment applies to a rotor 
that has any number of blades greater than tvro. Only a 
brief discussion and the formula for shaft-critical speed 
are given for the one- or two-blade rotor. 

The use of complex variables in conjunction with 
Lagrange's equations has been found very convenient for 
the treatment of vibrations of rotating systems. 



INTRODUCTION 



A rotary-wing aircraft that has hinged blades will, 
under certain conditions, be subject to vibrations which 



derive their energy'from the rotation of the rotor in- 
stead of from the air forces. The term "ground reso- 
nance" usually refers to vibrations of this type. Al- 
though such vibrations have apparently caused accidents 
in some rotary-ving aircraft and have impaired the flying 
qualities of others, very little attention has "been given 
this problem in the literature. A theoretical analysis 
has therefore been undertaken, and the purpose of the 
present paper is to present the theory and to describe 
the application of the theory to rotary-wing aircraft. 

General vibration theory and its application to al- 
lied problems as well as to the particular problem of 
rotor vibration are discussed in references 1 to 4. A 
good general background for the present problem is provid- 
ed in the chapters on rotating machinery and on self- 
excited vibrations in reference 1. Preferences 2 and 3 
treat in more abstract fashion the topics of rotation and 
damping. A discussion of the variety of modes of vibra- 
tion that exist in rotors and a number of frequency formu- 
las obtained by considering separately each degree of 
freedom are given in reference 4. This discussion does 
not, however, lead to a prediction of self-excited modes 
of vibration, 

Experience has shewn that two types of mechanical 
vibration may occur in rotors. The vibration frequency 
of the pylon is equal to the rotational speed in one type, 
unequal in the other. The first type is sometimes called 
the even-frequency vibration or the one-to-one frequency, 
and the second type, the odd frequency. The one-to-one 
frequency vibration resembles the phenomenon occurring at 
a critical speed of the shaft of rotating machinery and 
will consequently be referred to in this paper as a 
11 shaft-critical vibration." The odd-frequency vibration 
is properly called a self-excitad vibration. 

Jt derivation of the equations of motion for vibrations 
of a rotor for the case in which the pylon stiffness is 
equal in all directions of deflection is contained in ref- 
erence 5. The equation for the shaft-critical speed is 
obtained and checked by tests of simple models.. Reference 
5, however, contains incorrect statements regarding the 
existence of self-excited vibrations. The error was due 
to a confusion in the use of conjugate complex quantities 
which has now been cleared up. The present paper there- 
fore supersedes reference 5 and, in order t o make the 
present pap sr independent of reference- 5, the complete 



derivations are included herein without reference to the 
earlier report . 

An alternative derivation of the characteristic equa- 
tion for the whirling speeds of a three-blade rotor has 
"been given by Wagner of the Kellett Autogiro Corporation. 
By considering only the case of a pylon having equal stiff- 
ness in all directions of deflection, Wagner has shortened 
the analysis "by considering directly the equilibrium of 
forces and moments under conditions of steady circular 
whirling, Some examples of the dependence of whirling 
speed upon rotational speed are given, and the formula 
for the shaft-critical speed is obtained. 

In the present report, the theory is extended to in- 
clude the effects of damping in the hinges and in the huh 
and the effects of different stiffnesses of pylon deflec- 
tion in different directions. She method of analysis, 
particularly the use of complex variables in the equations 
of motion, is explained in some detail and all the previous 
results are shown to be a special case of the more general 
problem here treated. 



SYMBOLS 



1 1 



■11 



AA 



1 1 



^•31 

Ass 

A 22 



B 



radial position of vertical hinge 



•13 



A 



1 2 



= A 



31 



elements of the characteristic deter- 
minant (see equation (31 )) 



distance from vertical hinge to center of mass of 
blade 

damping force per unit velocity of pylon displacement 
/ E_ + B..\ 
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B x - B y 

£B = — 2- 

2 

Bj coefficient defined in equation (35) 

Br coefficient defined in equation (34) 

c,C lf ... C 4 arbitrary constants 

Cj coefficient defined in equation (36) 

C E coefficient defined in equation (34) 

D time-derivative operator (d/dt) 

F dissipation function 

I moment of inertia of "blade about hinge 

-b»"(i*S)] 

5 coefficients defined in equation (37) 



l » 



j s k indices and subscripts used with hinge coordinates 
(equation (l4)) 



K v + K, 



K spring constant (K f = J 



A. ^- — A.« 

Att _ t/ 
x , f = 



m 



/ m x + m,,\ 
effective mass of pylon ( m f = — 



m x - n y 
Am . _____ 

M total effective mass of "blades and pylon (m + nm^) 

AM mass added at hub for vibration test 
n total number of blades 
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r radius of gyration of "blade about its center of mass 

coefficients defined in equation (37) 
s stiffness ratio ( K y/^x) 

t time 

T kinetic energy 

T r kinetic energy of rotation of "blade about its center 

of ma s s 

T k kinetic energy of translat i onal motion of kth blade 
T_ kinetic energy of pylon 
V potential energy 

x,y displacements 

x 0 ,y 0 values of x and y when t = 0 

z complex displacement (x + iy) 

z complex conjugate of z (x-iy) 

a angle between blades (^jf) 

$ 0 , P x , ... P k angular displacements of blades 

P k value of P k when t = 0 

t 0 » ii* ••• £k variables representing hinge deflections 

when equations are expressed in fixed 
coordinate system 

G 0 , Q 1 , ... 6 k variables representing hinge deflections 

when equations are expressed in rotating 
coordinate system 

v == _x ( — x_ i n a p-piicat ions ) 

X M x \H x Wr J 

in applications ) 



My VMyWj 
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= ^j- (yj~ in applications^ 



*a 



Ax = 



A 



U) 

r 



M 

B _a 

M 



M 



0 

A a , = ( y^"2" in applications^ 



mass ratio 



, nm 
Ve + 



nm D y 

u i , u 2 expressions defined in equation (3) 

O) angular velocity of rotor (the dimen s i onle s s ratio 

uj/vJ r is called uj in applications) 

a> & angular whirling velocity measured in rotating co- 

ordinate system (used in nondimensional form in 
appli cat ion? ) 

iUf angular whirling velocity measured in fixed coordi- 

nate sysr,ets '■used in nondimensional form in ap- 
plications) 

reference frequency ( yK x /M x ) 
Suhscript s : 

f fixed coordinate system 

a rotating coordinate system 
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hinge deflection 



x,y component directions in fixed coordinate system 
b "blade 

APPROACH TO THE VIBRATION PROBLEM 
Stability and Instability 

If a vibrator were attached to a rotorcraft, several 
modes of vibration could be excited at any rotor speed. 
Only the modes that are likely to be excited during oper- 
ation of the aircraft, however, are important. 

In the present discussion, it is convenient to classi- 
fy the modes of vibration according to the circumstances 
required for their excitation. The different types of vi- 
bration are identified analytically by the nature of the 
roots of the characteristic equation. A hinged rotor may 
encounter three types of vibration which, for convenience, 
are herein designated ordinary, self excited, and shaft 
critical. At zero or slow rotational speeds, an external 
force is required to excite vibration. The friction al- 
ways present in such systems causes the vibration to be 
damped out when the force is removed. Modes of vibration 
requiring an external applied force to maintain them will 
be called ordinary vibrations. The mathematically ideal- 
ized case of zero damping will also be considered an ordi- 
nary vibration when it is understood to he an approxima- 
tion to a system actually damped. Self-excited modes of 
vibration are those with negative damping and are recog- 
nized analytically by the fact that a root of the charac- 
teristic equation is a complex number which has a negative 
imaginary part. A slight disturbance will tend to in- 
crease with time instead of damping out . 

When a rotating system is not perfectly balanced, the 
centrifugal force of the unbalanced mass may excite vibra- 
tions that have peak amplitudes at certain rotational 
speeds. Vibration excited by unbalance and in resonance 
with the rotation will be called shaft-critical vibration. 
This type occurs at the rotational speed at which the 
spring stiffness of the pylon is neutralized by the cen- 
trifugal force. In the analysis, the shaft-critical vi- 
bration is recognized in rotating coordinates as a zero 
frequency and in fixed coordinates as a frequency equal to 



8 



the rotational speed. The critical speeds of a rotating 
shaft are a common example of this class. 

The purpose of a theory of rotor vibration is mainly 
to predict the occurrence of and, if possible, to show how 
to avoid self-excited and shaft-critical vibrations. 

Choice of Degrees of Freedom 

Of the large number of degrees of freedom of a hinged 
rotor, the important ones for the present problem have 
been found to be hinge deflection of the blades in the 
plane of rotation and horizontal deflections of the pylon. 
Other degrees of freedom, such as the flapping hinge mo- 
tion of the blades, the bending or torsion of the blades, 
and the torsion of the drive shaft, are considered unim- 
portant in the problem of self-excited oscillations. Some 
motions, such as landing-gear deflection, that produce 
lateral deflection at the top of the pylon may, however, 
be imp or t ant . 

Physical Parameters 

She theoretical results given later provide a means 
of predicting the natural frequencies and, in particular, 
the critical speeds and unstable speed ranges in terms of 
certain physical parameters, such as mass, stiffness, and 
length. The successful application of the theory depends 
upon the determination of the proper values of these phys- 
ical parameters for the aircraft or model being studied. 

The important parameters t o be determined are: 

a radial position of vertical hinge. 

b distance from vertical hinge to center of mass of 

blade . 

m-jj mass of blade. Flexibility of the blade structure 
may have to be taken into account by tho use of 
an effective value of m-jj different from the ac- 
tual blade mass. The effective blade mass can be 
taken as the value required to make the theory 
predict the correct pylon natural frequency when 
the rotor has a. zero or very slow rotational speed. 

I moment of inertia of blade about hinge m^b 5 (l + h^) 
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00 



Kg spring constant of self-centering springs, which can he 

determined hy a force test or from the hinge frequency 
with the hub rigidly supported. 



& m^my effective mass of pylon for deflections in x- and y-directions. 

^ KzyKy effective stiffness of pylon. 

The effective mass of the pylon is the value of a con- 
centrated mass that would have the seme kinetic energy expressed 
in terms of the deflections at the hub as the actual pylon and huh 
if it were placed at the rotor hub in the plane of rotation. 
The effective ^stiffness of the pylon- is the stiffness of a 
spring that, if placed at the hub in the plane of rotation, would 
have the same potential energy in terms of deflections at the huh 
as the actual pylon. Equivalent definitions are that, if a simple 
spring and mass were imagined to be substituted at the hub in the 
plane of rotation for the pylon and hub, the natural frequency and 
the change of natural frequency with added mass would be the same as 
for the actual pylon. 

An experimental method of measuring the effective mass %- 
and stiffness K x of the pylon is to replace the rotor hy an 
approximately equal, rigid, concentrated mass £M at the huh 
and to measure the natural frequoncy for two or more values of this 
added mass. The quantities are then related hy the equation 



E 

03f = 



or 



0^2 K x 



If measured values of l/cof 2 are plotted against added mass AM 
and a straight line is drawn through the points, the intercept and 
the slope of the line will determine the effective values of 
Kjc and m^ . 

For the parameters a and h, the actual geometric lengths should 
he used unless the flexibility of the hinge off set arm a is 
comparahle in magnitude with the hinge spring stiffness. In this 
case, it is recommended that an effective value of a be guessed 
and that h be determined hy subtraction from the correct geometric 
value of a + b. 
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The damping parameters may "be defined "by the form of a 
dissipation function P. The function 

, p ,p S^l , 0 

2F a B z x f *= + Byyf + _> Bpfijj. 

K=0 

is equal to the rate of dissipation of energy hy damping. 

The parameters B x and By thus measure the damping force per 

unit velocity referred to linear displacements of the top of the 
pylon and Bp is the damping torque per unit angular velocity 

at a blade hinge. If the damping force per unit velocity is not *. 
constant, effective values should he used that will represent the 
same dissipation of energy per cycle as actually occurs with a 
reasonahle amplitude of vibration. The amplitude of free 
vibration in a single degree of freedom is given in terms of 
B x , B y , and B p hy 

Xf * Xq 0 ^ 

o 



The damping parameters are prohahly the most difficult ones 
to measure accurately. In practice, it is advisahle to make 
calculations for a given case, first on the oasis of no damping 
and then with the use of the estimated values of the damping 
parameters. 



MATHEMATICAL ISimOPMEHTS 
Method of Analysis 



The derivation of the characteristic equation that is used 
as the basis for predicting the unstahle oscillations of a rotor 
is presented in this section of the report. 
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Readers interested solely in. applications can omit this 
section and proceed immediately to the section entitled 
*i*fethod of Applying Theory." 

The method of analysis treats the equations of motion 
for small displacements from the equilibrium condition 
with steady rotation. A proper choice of coordinates 
leads to equations with constant coefficients. The solu- 
tions are exponential or trigonometric functions. 

The mathematical manipulations involved in treating 
the motions of a mass in a plane of rotation are facili- 
tated by the use of a complex variable. The typical dis- 
turbed motion obtai&sd by solving the equations of mo- 
tion is an elliptic whirling motion, which is represented 
in terms of a complex variable z = x 4 iy. At any in- 
stant, z represents the displacement of the pylon from 
its equilibrium position. An expression such as 



ce 



iou-f t 



represents whirling of the pylon in a circle of radius < 

with angular velocity u) f . The sign of ujf determines 

the sense of the rotation. Two rotations in opposite 

sense with the same radius are equivalent to a vibratory 
motion in a straight line. 



/ iuu-t 
z = c f e 1 + e 

= 2c cos oift 



-ioi^t 



f^ 



Two opposite rotations of different radii are equivalent ti 



whirling in an ellipse. A complex value of 



U) 



f represents 



whirling in a spiral, which may be either a damped or a 
self-excited motion depending upon the sign of the imagi- 
nary part. A self-excited motion exists when the imaginary 
part of u) f is negative, and the magnitude of z increases 
with time. 

The displacements may be referred to a fixed or to a 
rotating coordinate system. If Zf and z a are the dis- 
placements with respect to a fixed and to a rotating ref- 
erence system, respectively , 



Zf « z a e 



iu>t 



If 
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then 

8 f 

A whirling speed U) a with respect to the rotating coordi- 
nates thus corresponds to a whirling speed U)f = ou a + U) 

with respect to the fixed coordinates. A shaft-critical 
vibration corresponds to uu = 0 in the rotating coordi- 

nate system or to u) f = ud in the fixed coordinate system. 



= ce 



= ce 



i (uu a +a))t 



Example of E-otor with Locked Hinges 

An example that involves a partial use of complex 
variables is given on page 253 of reference 2. The prob- 
lem given there of a mass particle moving on the inner 
surface of a rotating spherical bowl is mathematically 
equivalent to the disturbed motion of a flywheel and shaft 
or of a rotor with locked hinges. The equations of motion 
obtained in real form in rotating coordinates 



S a - 2u# n - w 2 x a = - 



+ 2uux Q - 



a 



w 3 y a 



M 



M 



K 

M X£ 



H V& 



(1) 



are combined in the single equation 

(z\[X> + ^\ z 
V M / 



5 0 + ( 2iiu + z 0 + ( - - w*) z a =' 0 



(2) 



where 



z a = x a + i^a 

is the complex position vector in the rotating coordinate 
system. The complete solution, if small damping is as- 
sumed, is 



wher e 
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1 B a / W\ 



r- 

1 B. 



v„ = ( 1 + — ) 

3 2 M \ wj 



The path of the motion is represented by rotations of a 
complex vector in a plane. 

The use of a complex variable has thus cut in half 
the number of equations to he handled and has yielded a 
solution from which the geometric path of the motion may 
easily 'be reconstructed. The advantage of the z-notation 
is not fully realised, however, unless it is used from 
the very beginning of the problem. The close similarity 
of this problem to the rotor-vibration problem makes it 
worth while to show the full application of the z-notation 
to the preceding example. The complex variable e a at 
any instant determines the position of the mass particle 
relative to the rotating coordinate system. If the posi- 
tion in a fixed coordinate system is denoted by Zf, 

z f = Za e i( ^ (4) 

and z a can be treated as a generalized coordinate in the 
Lagrangian equat i on-s ; ~"'Of motion. The kinetic and potential 
energy expressions can be immediately written as 

1 „• •-*- 

= - M(z a + iujz a ) (z a - iujz a ) 
Y = i K Za z a 

A dissipation function for damping that depends upon motion 
relative to the rotating system can be written 

F " \ S a^a 

The equations of motion are now obtained by considering 
z a and ~z a as generalized coordinates in the Lagrangian 
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equations. Substitution in the equation 



\.-4 7 d2 a Bz a Bz a 



•thus yields the equation previously given 



"he same method can be used to obtain the equations 
of motion in the fixed coordinate system. In thir. case, 

"I 



V = 



- K« f T f 



The equation of motion in terms of Zf becomes 



Zf + 



Bo / • \ K 



and the solution for small values of damping is 



= C x e 



.. fl _ _J£„_ V '+ i /kTm t 



B 



+ C 3 e 



a 

2M 



\ (i + ^ ") t - i TkTm t 



(5) 



(6) 



(7) 



This solution shews that the motion consists of two circu- 
lar vibrations In opposite directions and, moreover, that 

for u) > </~T/U 'blio fi:-.-s-.i tc:--a :/e;p:. cscn,s '3'v:le action* blict is, 



This example illustrates a shaf t- cr i t i cal speed for 
uu = JlijM and a self-excited instability for uu > y K/M. 
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on 

i 



A discussion of the physical picture of this instability 
due to damping is given on page 293 of reference 1. 

The effect of damping in a nonrotating part of the 
system can "be included in the analysis merely "by adding to 
the previous dissipation function the term 

1 • _L 

2 B f z f z f 



The equation of. motion then "becomes 



Mz f + B 



\ 



f + B a (*f ' i{X)Z f) + Xz f = 0 
The solution for small values of damping "becomes 

'3.ifi. a ^ + i jtm] t 



(8) 



+ 0 2 o 



5* 

2M 



2M 



(i +-7=) - 



The motion is now unstable above the speed 



K 



O) = 




M (1 + 



B. 

b1 



(9) 



(10) 



Hinged Rotor 

Inclusion of the effect of hinge motion in the plane 
of rotation increases the number of degrees of freedom 
and the number of equations of motion. For example, three 
hinged blades and two directions of pylon deflection give 
five degrees of freedom to be considered. If special lin- 
ear combinations of the hinge deflections (3^ are used as 
generalized coordinates, no more than four degrees of free, 
dom need be considered simultaneously. The use of complex 
variables reduces these four equations to two equations. 

Appropriate variables in the rotating system for a 
three-blade rotor are 
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Go = 



M 
3" 



(p 0 + Pi + Pa) 



2-rri 



4TTi 



bi ! ° 3 



3 



e a = 



4TTi 8TTi 

/ 3 3 

iQ c + Pie + P 2 e 



(11) 



These variables and their complex conjugates satisfy the 
relations 



and also 



e x - - e a e a = - e x 
Mo + + = ~ e o 3 + ae^i - ^ ( p o 2 + P, 2 +■ P 3 3 ) 

The variables P k , by virtue of their meaning, are re- 
ferred to a rotating. coordinate system. The special linear 
combinations of the P^ denoted "by 6^ are also referred 
to a rotating coordinate system. The appropriate variables 
to represent the hinge deflections when fixed coordinates 
are used are defined by 



£k = e k e 



iuut 



(12) 



Geometrically, 8 X or £i is the complex vector rep- 
resenting the displacement due to hinge deflection of the 
center of mass of all the blades, just as z represents 
the position of the shaft due to pylon deflection. It will 
be shown later that, in the equations of motion, Q 1 is 

coupled with z and 8 Q is an independent principal co- 
ordinate. Equations (ll) when solved for 0 O , P lt and B 3 
become 



e 0 + 



+ e. 



= bip, 



6 0 + 6 x e ia + G 3 e ia = biP 1 



6 0 + G ie ia +6_e- ia = biP. 
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Then, in a mode involving 8 x i 





Qo = 


0 






6i = 


b e iw a t 

2 






e 3 = 


- e x 




^0 


= sin 






Pi 


= sin 


(u) a t - a) 


> (13) 


p 3 


= sin 


(a> t + a) 





Sanations (13) show that, in the Gi -node , the blades 
are undergoing sinusoidal vibrations 120° out of phase with 
one another in a manner analogous to three-phase electrical 
current s . 

General formulas for any number of blades are 

n-i 




(14) 



k=0 k=0 j 



Derivation of Equations of Motion 

The equations of motion and the characteristic equa- 
tion of whirling speeds are herein derived for the general 
case of three or more equal blades on a pylon that may 
have different stiffness properties in different direc- 
tions of deflection. The effects of damping in the blade 
hinges and in the pylon are included. The equations are 
first formulated in a nonrotating reference system. The 
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required modifications ere then given for the case of isotropic 
support stiffness. The corresponding equations referred to the 
rotating coordinates are theja obtained. 



Let the position of the center of mass of the kth 
"blade "be represented by the complex quantity z^. in the 
plane of rotation. (See fig. 1.) Let the "bending deflec- 
tion oi -he pylon he represented "by z^ in a nonrotating 
coordinate system and let 0 k he the hinge deflection of 
the kth blade. Then 



iPlc\ i(ak+uut) , . 

» k = z f + fa + be ) e (15) 



The complex velocity is 

z k = a £ + bi3 k e 1Pk + iu) (a + be iek ) 



i ( ak+cut ) / _ \ 
e (16) 



Because only small displacements are being considered, the 
exponential factors containing 3 k can be expanded and 

only the terms that lead to quadratic terms for the kinetic- 
energy expression need be considered. 

Some terms can be ignored either because they cancel 
after summation for all the blades or because the corre- 
sponding derivative expressions in the Lagrangian equations 
vanish. The substitution 

e * = 1 + iP k - -f- 

leads to an expression for the kinetic energy of transla- 
tional motion of the kth" blade. 

• T k = \ m b zk "2k (l?) 

where 
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The kinetic energy of rotation about the center of mass of 
the blade is 



-r = 3 *b r 2 p k * 



(18) 



The effective mass of the pylon may he different in the 
x f - and in the y f -directions . Allowance for this possi- 
bility is made by writing the kinetic energy of the pylon 
as 



wh e r e 



T _ = 



1 / • 3 

1 vVf 



= — j ms z _ + Am 



— i m s z 
2 [ f f 



( 



)] 



(19) 



m = 



m„ + m„ 
x J 



Am = 



m x ~ m y 



The total kinetic energy is the sum of the expressions for 
the separate kinetic energies. 



The pylon spring constant may differ in the x f - and in 
the y^-directi ons and, consequently, the potential energy 
can be ex-pressed as 



V = 



3,-2 

Z f + z f 

K Zf ? f + AK — + 



n- 1 



p p k 



k=0 



(20) 



The effect of damping will be expressed with the aid 
of a dissipation function. If damping exists in the pylon, 
in the rotating shaft, and in the hinges, this function be- 
c one s 



F 



i 



f z f 



„ 2 . — 2 
f f ' ' 

+ AB . _ + B aZa z a + 



n-i 



Vk 



k=0 



(21) 



The sum of the various energy expressions for all the 
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blades, expressed in terms of the variables z f and £ k 
in the nonrotating coordinates., becomes 



Z Z. 2 

. + Z~ • J_ 

' + (n+nn^) z f z f + 



nQ b ( z f L + »fL) 



1 

V = — 

2 



• 2 + 7. 



AS 



> (22) 



1 

2 l 



z^. 2 + z f 2 • • 

A3 — =— + B Zf z f 



+ B a ^ Zf-io)Zf^ z f +ituz f ^ 

^ y (t k - ^0 (i k - i-to] 



n-l 



b ^ 

The La^rangian equations of motion are 



J 



_d 

dt ^31 



= o 



3 Zf 3 z f 3 Zf 



1_ ,'3T _\ 3T_ + aP_ + cTV_ = 



(23) 



dt Vo^/ 



and similar expressions for the other variables. The equa- 
tions of motion in fixed coordinates then become 

(m+nmfc) i f + Bz f + 3 a (z f -iUiz f ) + Kz f + Amz" f + ABt f + AEz f + nn^ ^= 0 



<2l+) 



-3 s k 



=r 0 
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where £ k refers to the ^-variables other than £1 . The 
complex conjugates of these equations are also obtained 
hut give no additional information. Each complex equation 
is, of course, equivalent to two real equations. It is 
°° noticed that the first two equations contain only the var- 

1 iables z_, "z_, and £ x and that the third equation 

represents n~2 equations, each containing one independ- 
ent principal coordinate £ k . The physical meaning of 

this partial separation of variables is that a blade mo- 
tion represented by £ x involves a motion of the common 
center of mass of the blades and, thus, a coupling effect 
with the pylon. Blade motions in which the common center 
of mass does not move are represented by £ 2 , ... £ n . For 
three blades, the only such mode is the one corresponding 
to £ 0 . In this mode, all the blades move in phase; the 

motion is always damped and does not lead to instability. 

The equations of motion of a one- or two-blade rotor 
are somewhat different from equations (24). The differ- 
ence is connected with the circumstance that a rotor of 
three or more equal blades has no preferred direction in 
its plane; whereas, a one- or two-blade rotor has differ- 
ent dynamic properties in directions along and normal to 
the blades. Only a brief statement- and the final equation 
for shaft-critical speed will be given for the one- or 
two-blade rotor. 

The equations of motion involving Zf and £j can 
be written more compactly by use of the notation 

D = — D 2 = d * 



dt at- 



and the substitutions 



B f Bp 

T = Af — ~~' ^ = A P 

mvb 



+ b~0 



B a . , 



M~ s A a r^~\ " A i 

* C 1 + jO 
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= A. 



!hen 



|~D a + * f D + A a (D-iuj) + | "L + D a + AX f L + -f ) 7 f + HD* £ x = 0 



> (25^ 



2 (l + O 



IT z f * 



(D-iw) 2 + >vp(D-iuj) + u> a A x + A £ 



= 0 



(26) 



or, "briefly, 

An (Ifcf + AA u Zf + A 1S (DXi = 0 
A 31 ^D) Zf + A 33 (r| 1 = 0 

The Characteristic equation 

The general form of solution of equations (26) is an 
elliptic whirling motion that can "be represented by 



J 







+ C^e- 1 ^ 






¥ f 


= "C 1 e" i '^f * 


+ C e 1 ^* 




> (27) 




= C 3 e iaJft 


iujf t 
+ 0 4 e 1 


— < 





Special cases of this motion include whirling in a circle 
C 2 = G 4 = 0, and linear vibration- c i = c s» C 3 = C 4 • 
Substitution of equations (2?) in equation (26) gives 



A xl (iuj f ) C 1 + AA X1 (i Wf ) C 3 +A 12 (iu) f ) C 2 
+ j^A lx (-iUJf) c z + AAnC-iuDC! + A 13 (-iu) f )C 4 



.ioif t 



e -iu) f t = o 



,io) f t 



j^A 13 (icu f ) Cj + A S2 (ioj f ) C. 

+ j^A 12 (-iuT f ) 0 2 + A aa (-iu»f) C 4 



> (28) 



e = 0 
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In order for equations (2?) to "be a solution of equations 
(26), equations (28) must "be satisfied for each value of 

ill) ft -iu3ft 
t. The coefficient of each time factor e or e 

must therefore separately vanish. Because each bracketed 
expression represents a complex quantity that vanishes, 
its complex conjugate also must vanish. The condition 
for a solution can therefore he expressed hy the vanishing 
of the first bracketed terms and the complex conjugates 
of the second bracketed terms. Hence, 



i u (iu) f )Ci +AA. 11 ( iou f )C 3 +A 13 (iuj f )C 2 



A 12 (i UJf )C 1 

AA 11 (i U)f )C 1 +A 11 (iuJ f )C a 



rA 
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(iu) f )Cj 



•is 



( iu)f ) 0. 



= 0 

= 0 

+A ia ( iujf )C, 4 = 0 

+I 32 (io; f )C 4 = 0 



} (29) 



.J 



whore A lx (iujf) is the complex conjugate of A 11 (-iuj ; f) 
and is obtained from A 11 (ioj f ) "by changing iw to -icu 
without changing iiiUf. The characteristic equation giving 
the rotational speeds is the determinant of the coefficients 



of 



i > 



0 S> C 3 , and C 4 



equated to zero. With the second 



and third columns interchanged for symmetry, the determi- 
nant becomes 





An 


A 1S 


AA 1X 


0 












0 


0 


= 0 


(30) 




AA 13 _ 


0 


X n 


A 1 








0 


0 


A i 3 


A 22 






The expanded form 


of this 


determinant 


is 






( A l 1 A 22 ~ A 12 A 21 ) 


^ A ll A 22- 


W 21 > 




&As = 0 


(3D 



where 
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A u = -aif 2 + iuufAf + iX a (aif-uj) + — 

A X1 = -a> f 2 + iu> f A f + iX a (uj f +uj) + | 

AA 1X = AA^ = - $ a> f s + ilWf AA f + 

A ia A 3 i = = — u> f 4 = A 3 w f 4 

2 C 1 + to 

A 33 = - (w f -cu) + iXp(o) f -o)) + (ju 2 ?^ + A 3 
A 22 = - (ujf+uj)"* + iXp(o> f +(Ju) + . , w a A 1 + A 2 

The roots ui^ of this equation are tha characteristic 
whirling speeds of the rotor. 

For the case of isotropic supports, 

aa i:l = 0 

*r>il the equations of notion are satisfied "by equations (27) 
with C 2 = C 4 = 0. 

The characteristic equation is then simply 



A n A 22 " Wsi = 0 (32) 

In a rotating coordinate system, the complex coordinates 
are z„ and 6, i where 

ieut 

z f = z a e 

Then 

iujt 



D Zf = (Dz a + iuiz a ) e' 
Dt x = (D6 1 + ioi6 x ) e ia,t 
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If the whirling speed in rotating coordinates is represent- 
ed by 0J a , 

z a = °i e a 

Q icu a t 
6 i - °s e a 

The characteristic equation is then obtained "by substitut- 
ing uj a + uj for uj f . 

A u (u) a +oi) A 32 (u) a +uu) - A 12 (uj a +u>) A ai ((ju a +aj) =0 (33) 

The characteristic equation can thus he stated in terms of 
a whirling speed in either the fixed or the rotating co- 
ordinate system., 

. METHOD OF APPLYING THEORY 
Application Neglecting Damping 

In plotting curves for use in applications of the the- 
ory, it is convenient to consider one of the pylon "bending 

frequencies o> r = J K x /M x as a reference frequency and 

to refer all other frequencies as well as the rotational 
speed uj to the reference frequency as unit. The number 
of independent parameters is thus redticed by 1. All 
quantities in equations (3l) to (33) are then expressed 
nondi men si onally . 

The natural whirling speeds and the three types of 
vibration - ordinary, self-excited, and shaft- crit i cal - 
can now be predicted from a study of the roots of equation 
(31) in which u> f is considered a function of U) for 
fixed values of the other parameters. 

The case of no damping will be considered first, Be—- 
cause equation (31) with damping terms omitted is of the 

2 

fourth degree in u) f and of only the second degree in 

uj 2 , it may be solved conveniently by first choosing val- 
ues of u)f and then solving the equation for u> 3 . Simi- 
lar indirect methods can be used with equations (32) and 
(33). Special methods to be used when damping is included 
will be discussed later. 
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The meaning of equations (3l) to (33) will "be illus- 
trated by examples. The real part of U) f will he plotted 
against u> for selected values of the parameters A x A 3 . 
A 3 , and e „ The simplest case is that in which the mass 

of the "blades is so small that any force on the pylon due 
to blade motions is negligible. The pylon motions are 
then independent of the blade motions. This case is ob- 
tained by putting As = °« iJ!ile characteristic equation 
(31), (33), or (33) then factors into expressions yielding 
straight lines and hyperbolas. 

An example of a rotor with particular values of the 
parameters is plotted as long-dash lines in figure 2. 
The horizontal straight lines correspond to pylon bonding 
and the slanting hyperbolas corrrespond to hinge deflec- 
tion. Bach curve represents the trend of one of .the real 
roots uv. As A 3 incroases slightly from zero, the 
greatest changes in the curves occur in the vicinity of the 
intersections of the straight lines with the hyperbolas. 
Here each branch breaks away from the intersection and re- 
joins the; other branch. At a gap, such as C in figure 
2, the number of real roots of the frequency equation is 
reduced by 2. The missing roots arc ccm.pl.ex conjugate 
numbers; and one of them must have a negat.ive imaginary 
part, which implies a seli-excited vibration. 

Consider the interpretation of figure 2 as a' is 
.gradually increased from zero. At zero rotational speed, 
the values of oi^ are the natural frequencies that could 

be excited as. ordinary vibration by applied vibrating 
force. Positive and negative values occur in pairs of 
equal magnitude and correspond to linear, vibration modes 
represented in complex notation as 



= c 



/ iu)f t - iuj-p t\ 
(o ; ♦ c, ) 



As uj increases from zero, the positive and negative 
valv.es of uj^, no longer are equal in magnitude. The nor- 
mal modes are therefore whirling motions with angular ve- 
locities ■■ equal t o the plotted values of oi£ . 

The shaft-critical speed is the rotational speed at 
which = u> and hence , is given by the point A where 

a 45° line through the origin intersects the ui^-curve. 
This speed corresponds to the peak for vibrations excited 
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"by unbalance in the rotating system. As u) increases 
above the shaft-critical speed, the modes of whirling are 
stable until, for the case of no damping, the value of o>£ 
"becomes complex at the value of u) at which a vertical 
line is tangent to the plotted curve. This point B is 
the "beginning of the self-excited range. At the point D, 
the motion again "becomes stable. She real part of UL' f is 

plotted in the region 0 as a short-dash line, The com- 
plex roots in the region 0 have been calculated and 
plotted in figure 3. 

The point 1, at which uu f = 0, is of some interest 

At this speed, a vibration of the blades could be excited 
by a steady force (u; f = 0, uj & = -u>) such as the force of 
gravity if the plane of the rotor is not horizontal. 

Because the most important information to be obtained 
from the frequency equation is the critical value of uu 
for the shaft-critical and self-excited vibrations, a set 
of charts that gives this information for a large variety 
of values of the physical parameters has been prepared. 
These charts are given in figures 4 to 6, which correspond 
to values of stiffness ratio K y /K x = s of 1 , « , and 0, 
respectively. The use of the charts is illustrated by a 
numerical example. Suppose the values of the parameters 
for a certain rotor are A x = 0.07, A s = 0.22, A 3 = 0.1, 

s = l, and u.' = 155 cycles per minute. A straight line, 
such as AB in figure 4, is first drawn to represent the 
function uu 3 A 1 + A a « This line intersects contours 
A 3 a 0.1 at a) 3 = 0.77 for the shaft-critical point and 

U) 3 = 1.6 and 4.85 for the beginning and for the end of 
the self-excited range, respectively. With a reference 
frequency of 155 cycles per minute, these values corre- 
spond to actual rotational speeds of 135, 195, and 342 rpm 

All possible values of A x , A 2 , and A 3 are thus cov- 
ered by suitably changing the straight line AB . The gen- 
eral effect of the stiffness ratio s is not large; any 
case can therefore be estimated with a fair degree of ac- 
curacy by use of figures 4 to S. 
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Possibility of Avoiding Occurrence of Titration 

Figures 4 to 6 can also be used for the inverse prob- 
lem of finding the values of the parameters that are re- 
quired to obtain given values of critical rotational speed 
These figures show that to eliminate entirely the self- 
excited instability requires that A x be equal to or 
greater than 1. The shaft-critical instability can be en- 
tirely eliminated only with a value of Ai in a small 
range near 4 and with s = a> or s = 0. These values of 
Ai differ radically from present designs in which a typ- 
ical value is 0.07. 

The satisfactory requirement of keeping the instabil- 
ities outside the operating range of rotational speed is 
found by first picking a reasonable value of the pylon 

frequency Jlt-fto to fix the scale unit for u) and by 

then obsarving the combinations of Ai and A a that can 

be used to avoid the critical A 3 -contours. 

Effect of Damping 

The effect of damping has been included in equation 
(31) through the parameters X x , A v , and Xfj . A method of 

computation similar to that used in flutter theory appears 
preferable to attempting to solve the equation directly 
for 0) fn The beginning and the end of an unstable range 

can bo found by the following method: At a limit point 
between a stable and an unstable speed range, the value of 
U) f is real. Equation (3l) is first separated into real 

and imaginary parts with iu^. considered real. Each part 

is considered a functional relation between ujf and uu 
and is plotted for a given set of values of the parameters 
The intersections of the real and the imaginary equations 
give the rotor speeds and frequencies corresponding to the 
beginning and the end of the unstable ranges. In the com- 
putations, it is preferable to choose values of uuf and 
to solve the equations for the corresponding values of. u) . 

The explicit form for computation in the simplest case 
isotropic supports.) wit'.; damping in the pylon and in the 
hinges but not in the rotating shaft (X a = 0 ) , is ob- 
tained from equation (32) rearranged as follows: 
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o 

"I 



Real equation 

U) 2 - 2Beuj + Or = 0 



where 



B* = 



1 - A, 



1 + 



2 (V + 1) 



and 



'E 



<JL) f « 



1 - A 3 



M X 



_U) f 3 + 



where 



and 



(34) 



(35) 



Imaginary equation 

uu 2 - 2BjU> + Cj = 0 

•"■•rhrKK-S^K.*'.)] 

The most general case obtained from equation (31) is writ- 
ten: 



Real equation 
Coefficient of oi 6 

UJ 2 

1 



T^aU-Ai) 3 

Sid-Aif + A a 8 R 3 

R i S 3 - + x a E s ~ R 5 

R X H 3 - I l I a - R 4 + A 3 a iu f e 



(36) 
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Imaginary equation 
Coefficient of co 4 

CD 3 

1 



1^1 - A,) 2 +X a 2 I 3 



R 1 I 3 + E 3 I 1 + X a " 

E l I 2 + Vx ' X « 



(37) 



where 



R. 



a fa ♦ SV- 



+ 3Ez ) -a 



a ilV M 



R 3 ■ ~ 2 ( X - A x) + A 2 ) " ^ + V 

E 4 " 4A 3 (^f 2 + |) ("^f 2 +A a ) "™f 2 (*f+*a) Xp j 
I 4 = ao f 4 A 3 j"^ f 2 + g jc»fXp + oi f ^X f + \ a ) (-o) f 2 + A^ 
r 5 = a, f 4 A 3 [- (l - Ai ) (-cn/ + |) -x a x p ] 

I s = 9o) f 4 A 3 £■ (l -A J oif (x f + x a ) -2X^1 

Examples of calculated cases with damping are shown in 
figures 7 to 9- The presence of small amounts of damping 
in both the pylon and the hinge degrees of freedom does 
not greatly change the predictions that would he made 
from the equations with no damping. The plot of the real 
equation is practically the same as the plot obtained when 
damping is neglected. The intersections of the curves of 
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the imaginary and the real equations with any reasonable 
value of ^f/^p are near the points that would he con- 
sidered the limits of the unstable range if damping were 
neglected. Increasing the amount of damping decreases 
the gap between the limits of stability until the unstable 
range is finally eliminated. An approximate solution for 
the amount of damping required to eliminate the self- 
excited instability is obtained by requiring that the 
damping be at least large enough to make the curve of the 
real equation pass through the point where <JUj = 1 and 

oi is the value given by the equation 



1 = a) - J u> 2 A x + A a 

The values required in the case of s = » have "been com- 
puted and plotted in figure 10. The elimination of self- 
excited vibration by damping thus looks promising and 
merits further study with reference to specific applica- 
tion. 

LIMITATIONS AND FURTHER DEVELOPMENTS OF THE THEORY 

Polar Symmetry 



An important idea in the rotor vibration theory is 
the concept of polar symmetry. This concept implies the 
absence of a preferred direction in the plane of the ro- 
tor. A rotor of three or'more equal "blades has polar 
symmetry. A rotor of two blades or one with unequal cen- 
tering springs does not have polar symmetry. A pylon 
for which K x = K y , B x = B y , and m x = m y has polar 

symmetry. The possibility of solving the rotor vibration 
problem in terms of exponential or trigonometric func- 
tions depends upon the existence of polar symmetry in the 
rotating parts or in the nonrotating parts or in both. 
The general case of no polar symmetry would lead to 
Mathieu functions or something similar. 



Two Blades 

A brief comparison between the two-blade and the gen- 
eral case is presented herein. Polar symmetry of the py- 
lon is assumed. The shaft-critical speed is obtained by 
substituting u). = 0 in the characteristic equation as 
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expressed in a rotating coordinate system. For one? or two 
blades, the equation obtained is 

f-V + ^ (fu; 3 ± + - H 0< 4 ] f-^ + = 0 (33) 

The first bracketed factor gives the beginning of a self- 
excited range and the second factor gives the end of the 
range . 

Equation (38) can be compared with the equation for 
the shaft-critical speed of three or more equal blades 
and for polar symmetry 



(-U, 2 + -) U - + ~\) - * ^ - 0 (39) 
V M/ V b m v b 3 / 2 



A useful chart based on equation (39) is given in figure 
11; some experimental results of tests of a simple model in 
figure 12. These tests demonstrate the essential differ- 
ence between the two-blade and the general case. 

Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, 7a., 
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Pigure 2.- The effect of coupling between pylon and hinge motions. 
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Figure 3.- The complex frequency in the unstable range for = .07; 
Ag = .22; A3 = .1; s = 1. 
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figure 4.- Stability chart for s = 1. 
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figure 5.- Stability chart for s = oo . 
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